Prime numbers with 
Beatty sequences 

William D. Banks 
Department of Mathematics 

University of Missouri 
Columbia, MO 65211 USA 
bbanksSmath.missouri . edu 

Igor E. Shparlinski 
Department of Computing 
Macquarie University 
Sydney, NSW 2109, Australia 
igorOics .mq. edu. au 

February 5, 2008 



Abstract 

A study of certain Hamiltonian systems has lead Y. Long to conjec- 
ture the existence of infinitely many primes of the form p = 2 \_an\ +1, 
where 1 < a < 2 is a fixed irrational number. An argument of 
P. Ribenboim coupled with classical results about the distribution 
of fractional parts of irrational multiples of primes in an arithmetic 
progression immediately imply that this conjecture holds in a much 
more precise asymptotic form. Motivated by this observation, we give 
an asymptotic formula for the number of primes p = q [an + /3J + a 
with n ^ N, where a, f3 are real numbers such that a is positive and 
irrational of finite type (which is true for almost all a) and a, q are 
integers with ^ a < q ^ N'^ and gcd(a, q) = 1, where k > depends 
only on a. We also prove a similar result for primes p = Ian + f3\ 
such that p = a (mod q). 
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1 Introduction 



For two fixed real numbers a and /?, the corresponding non-homogeneous 
Beatty sequence is the sequence of integers defined by 

Ba,i3 = (La^ + /5J)~=i- 

Beatty sequences appear in a variety of apparently unrelated mathematical 
settings, and because of their versatility, the arithmetic properties of these 
sequences have been extensively explored in the literature; see, for example, 
[D S H El m [m [121 [Tg [Tg [21] and the references contained therein. 

In 2000, while investigating the Maslov-type index theory for Hamiltonian 
systems. Long [T7] made the following conjecture: 

Conjecture. For every irrational number 1 < a < 2, there are 
infinitely many prime numbers of the form p = 2 \_an\ + 1 for 
some n eN. 

Jia [9] has given a lower bound for the number of such primes p in the interval 
{x/2,x]. We remark that, using a simple modification to an argument given 
by Ribenboim [20l Chapter 4.V], one can show further that the number of 
such primes p ^ x is asymptotic to a~^7i{x) as x — >• oo; see also [16]. More- 
over, Ribenboim's method also applies to the general problem of estimating 

N'a,i3-q,a{x) = ^ X : p = Q [au + P\ + a is prime}, 

where a, (3 are fixed real numbers such that a is positive and irrational, and 
a,q are integers with ^ a < q and gcd(a, g) = 1. In fact, if a and q are 
fixed, one easily derives the asymptotic formula 

q 

A/'a,/3;g,a(a;) = (1 + o(l)) ^^7r(x) (x^oo), 

where the function implied by o(-) depends on a, and q, and (p{-) is the 
Euler function. Motivated by this observation, we consider here the problem 
of finding uniform estimates for Afa,i3;q,a{x) if q is allowed to grow with x. We 
also consider the same problem for the counting function 

M.a,i3;q,a{x) = ^ X : p = [an + P\ is prime, and p = a (mod g)}. 

In particular, in the case that a is of finite type (which is true for almost all a 
in sense of Lebesgue measure), our main results yield (by partial summation) 
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nontrivial results for both J^a,/3;g,a{x) and ^Aa,i3;g,a{x) even if q grows as a 
certain power of x. 
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boim [20] on primes in a Beatty sequence, and Ahmet Giiloglu for pointing 
out a mistake in the original version of the manuscript. During the prepara- 
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2 Notation 

The notation is used to denote the distance from the real number x to 
the nearest integer; that is, 

||x|| = min \x — n\ (xE R). 

nez 

As usual, we denote by [x\ , \x] , and {x} the greatest integer ^ x, the least 
integer ^ x, and the fractional part of x, respectively. 

We also put e(x) = e^'^*^ for all real numbers x and use A(-) to denote 
the von Mangoldt function: 




logp if n is a power of a prime p; 
otherwise. 



Throughout the paper, the implied constants in symbols O, <^ and 3> 
may depend on the parameters a and P but are absolute unless indicated 
otherwise. We recall that the notations A <^ B, B ^ A and A = 0{B) are 
all equivalent to the statement that |y4| ^ c\B\ for some constant c > 0. 
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3 Preliminaries 



Recall that the discrepancy D{M) of a sequence of (not necessarily distinct) 
real numbers ai, a2, . . . , G [0, 1) is defined by 

V{I, M) 



D{M) = sup 

XC[0,1) 



M 



\I\ 



(1) 



where the supremum is taken all subintervals X = (c, d) of the interval [0, 1), 
V(X, M) is the number of positive integers m ^ M such that G X, and 
|X| = d — c is the length of X. 

For an irrational number 7, we define its type r by the relation 



r = sup|f) G R : limijtif ||7^|| = o| 



(nGN) 

Using Dirichlet's approximation theorem, it is easily seen that r ^ 1 for 
every irrational number 7. The celebrated theorems of Khinchin |TD] and 
of Roth [2TI [22] assert that r = 1 for almost all real (in the sense of the 
Lebesgue measure) and all irrational algebraic numbers 7, respectively; see 
also [SI [23]. 

For every irrational number 7, it is well known that the sequence of 
fractional parts {7}, {27}, {87}, . . . , is uniformly distributed modulo 1 (for 
instance, see [14| Example 2.1, Chapter 1]). If 7 is of finite type, this state- 
ment can be made more precise. Let D^^i{M) denote the discrepancy of the 
sequence of fractional parts {{'jm + S})^^^. By [Tl[ Theorem 3.2, Chapter 2] 
we have: 

Lemma 3.1. Let ^ be a fixed irrational number of finite type r < 00. Then, 
for a// (5 G M the following bound holds: 

D^^M) ^ M-i/"+°(^) (M 00), 

where the function implied by o(-) depends only on 7. 

The following elementary result characterizes the set of numbers that 
occur in a Beatty sequence Ba,f3 in the case that a > 1: 

Lemma 3.2. Let G M with a > 1. Then, an integer m has the form 
m = \_an + j3\ for some integer n if and only if 

< {a~^{m- (3 + 1)] ^ a~^. 

The value of n is determined uniquely by m. 
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Proof. It is easy to see that an integer m has the form m = \_an + (3\ for 
some integer n if and only if the inequahties 

m — 13 m — 3 + 1 

^ n < 

a a 

hold, and since a > 1 the value of n is determined uniquely. □ 

We also need the following statement, which is a simplified and weakened 
version of a theorem of Balog and Perelli [2] (see also [15] ) : 

Lemma 3.3. For an arbitrary real number and coprime integers a, q with 
^ a < q, if \d — h/d\ ^ 1/L and gcd(6, d) = 1, then the bound 

J2 Mn) ei^n) « (-^ + d'^'l'^' + L'A (log L)' 

n=a (mod q) 

holds, where the implied constant is absolute. 

Finally, we use the Siegel-Walfisz theorem (see, for example, the book [8] 
by Huxley), which asserts that for any fixed constant B > and uniformly 
for integers L ^ 3 and ^ a < q ^ (logL)^ with gcd(a, q) = 1, one has 

J2 A(ri) = -^ + 0(Lexp(-CBV^)), (2) 

n=a (mod q) 

where Cb > is an absolute constant that depends only on B. 



4 Bounds on exponential sums 

The following result may be well known but does not seem to be recorded in 
the literature. Thus, we present it here with a complete proof. 

Theorem 4.1. Let 7 be a fixed irrational number of finite type r < 00. Then, 
for every real number < e < l/(8r), there is a number t] > such that the 
bound 

A(gm + a) e{jkm) 

holds for all integers 1 ^ k ^ and ^ a < q ^ M^/^ with gcd(a, q) = 1 
provided that M is sufficiently large. 
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Proof. Fix a constant g such that 



1 ^ ^ < ^ < ^ 
8e 



(3) 



Since 7 is of type r, for some constant c > we have 

> cd-^ {d ^ 1). 



(4) 

Let k, a, q be integers with the properties stated in the proposition, and 
write 



A{qm + a) ei^jkm) = e{—{}a) A{n)e{'dn), 



(5) 



n=a (mod q) 



where = jk/q and L = qM + a. Let b/d be the convergent in the continued 
fraction expansion of ^ which has the largest denominator d not exceeding 
L^-^; then, 



7A; b 
q d 

Multiplying by qd, we get from (jl]): 



1 



dL 



l-e 



(6) 



L 



^ ^ \-fkd-bq\ ^ \\-fkd\\ > c{kd) 



Thus, since we see that under the condition ([3]) 

the bound 

d ^ CL^^-^'^/e-' ^ CL^'^^^^ (7) 

holds, where C = c^^^ and L is sufficiently large. 

Inserting ^ into ([6]) and using ([3]) again, we conclude that 

7fc 6 1 1 

if L is sufficiently large. We are therefore in a position to apply Lemma 13. 3|: 
taking into account ([3]), 1 
bound 



and the fact that d ^ L ^, it follows that the 



A{n) e{^n) < [L^-^i^^) + L^-''^) {\ogLf ^ L^' 



■e/3 



n=a (mod q) 



holds for all sufficiently large L. Since L <C qM ^ M^"*"^/"^, the result now 
follows from simple calculations after inserting this estimate into □ 
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Using similar arguments, we have: 



Theorem 4.2. Let 7 be a fixed irrational number of finite type r < 00. Then, 
for every real number < £ < l/(8r), there is a number 77 > such that the 
bound 



A(m) e{'ykm] 



m=a (mod q) 



holds for all integers 1 ^ k ^ and ^ a < q ^ M^^^ with gcd(a, q) 
provided that M is sufficiently large. 



5 Main Results 

Theorem 5.1. Let a and (3 be a fixed real numbers with a positive, irrational, 
and of finite type. Then there is a positive constant k > such that for all 
integers ^ a < q ^ N'^ with gcd(a, q) = 1, we have 



A{q [an + p\ + 



n<N 



a) = a-^ J2 ^(^"^ + a) + {N^-^) 

m^[aN+(3\ 



where the implied constant depends only on a and (3. 

Proof. Suppose first that a > 1. It is obvious that if a is of finite type, then 
so is a~^. We choose 

1 

< e < , 

16r 

where 1 ^ r < 00 is the type of a^^ . 

First, let us suppose that a > 1. Put 7 = 5 = a~^{l — P), and 
M = [aN + p\. By Lemma ESI it follows that 

Sa,l3;q,a 

(N) = A(g [an + p\+a) 

Y A{qm + a) + 0(1) 

0<{7m+5}<7 

Y + «) 'ipiim + 5) + 0(1), 



ri<Af 



(8) 
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where ip{x) is the periodic function with period one for which 

M = I ^ if < a: ^ 7; 
' \ if 7 < X ^ 1. 

By a classical result of Vinogradov (see [2Sl Chapter I, Lemma 12]) it is 
known that for any A such that 

1 1 
0<A<- and A ^ - min{7, 1 - 7}, 
8 2 

there is a real- valued function with the following properties: 

(i) i^Aix) is periodic with period one; 
(ii) ^ V^A(a;) ^ 1 for all X G M; 

{iii) ipAix) = ipix) ifA^a;^7 — Aorif7 + A^a;^l — A; 
(iv) iPa{x) can be represented Fourier series 



k=l 

where the coefficients satisfy the uniform bound 

max{|5(fc|, < min{A;"\ A;"^A"^} (k^l). (9) 

Therefore, from ([8]) we deduce that 

Sa,(3;UN)= Yl Hq^ + a)M7^ + S) + 0{l + V{I,M)\ogN), (10) 

where V{X, M) denotes the number of positive integers m ^ M such that 

{7m + 5} G J = [0, A) U (7 - A, 7 + A) U (1 - A, 1). 
Since |X| ^ A, it follows from the definition ([T]) and Lemma [3. Ij that 

F(T,M) < AiV + iV^^", (11) 
where the implied constant depends only on a. 
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To estimate the sum in fllUp . we use the Fourier expansion for ip^l'jm + S) 
and change the order of summation, obtaining 



A{qm + a) ^p/\{'jm + 5) 

oo 

7 A(gm + a) + e{6k) A(gm + a) e(7/cr 

fc=l m^Af 
oo 

+ /ifc e(— A(gm + a) e{—'ykm). 



zm] 



'12) 



m<M 



k=l 



By Theorem 14.11 and the bound Q), we see that for ^ a < g ^ M*^/^, 
we have 

^ e{5k) + «) ei'ykm) < M^"" ^ A;"^ < M^-^/^^ (13) 

for some ?/ > that depends only on a. Similarly, 

hke{-6k) Y Kqm + a) e{--fkm) < M^"''/^ 

On the other hand, using the trivial bound 



(14) 



Y^ A(gm + a) e{'ykm) 



n<N 



we have 



Y 9k e{Sk) J2 Ml^ + a) e{jkm) < iV ^ k'^A'^ < N^-'A-\ (15) 



and 



k>M^ 



Y hk e{-6k) Y Mq^ + a) e{--fkm) < A^^^^A^^ (16) 
Inserting the bounds, ([13]), ([IID, ([15]) and ([IS]) into we obtain 



A(qm + a) ^pAiim + S) 



=7 ^ A(gm + a) + 0{M^~'^l^ + A^^-^A"^), 



(17) 



m<A// 



where the constant imphed by O(-) depends only on a and (3. 

Substituting (fTTi) and (fT7|) in (fTUi) and choosing A = N~^^'^, it follows 
that 

^a,/3;,,a(iV) = 7 A(g^ + a) + O {N'-^) , (18) 

for some n which depends only on a. This concludes the proof in the case 
that a > 1. 

If a < 1, we put t = and write 

t-i 

^ A(g [an + + a) = ^ ^ A{q [atn + aj + P\ + a). 

n<Af j=0 n<i{N-j)/t 

Applying the preceding argument with the irrational number at > 1, we 
conclude the proof. □ 

In particular, using the Siegel-Walfisz theorem ([2]) to estimate the sum 
in (fT8|) for "small" a and q, we obtain: 

Corollary 5.2. Under the conditions of Theorem \5.1[ for any constant B > 
and uniformly for all integers N ^ 3 and ^ a < q ^ (logN)^ with 
gcd(a, q) = 1, we have 

^ A(g [an + /3J + a) = -p- iV + O (^N exp (^-C^/\ogN^^ 

for some constant C > that depends only on a, (3 and B. 

In the special case that (a, q) = (0, 1) or (1, 2) (the latter case corre- 
sponding to primes in the Long conjecture), we can use a well known bound 
on the error term in the Prime Number Theorem (proved independently by 
Korobov [I3] and Vinogradov [25]) to achieve the following sharper result: 

Corollary 5.3. Suppose that {a,q) = (0,1) or {a,q) = (1,2). Then, under 
the conditions of Theorem \5.1{ for any constant B > and uniformly for all 
integers N ^ 3, we have 

^ A{q [an + (3} + a) = q N + O {N exp (-c(log iV)^/^(log log A^)~^/^)) 
for some absolute constant c > 0. 
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Finally, using Lemma 14.21 in place of Lemma 14.11 we obtain the following 
analogues of Theorem 15.11 and its two corollaries: 

Theorem 5.4. Let a and (3 he a fixed real numbers with a positive, irrational, 
and of finite type. Then there is a positive constant k > such that for all 
integers ^ a < q ^ N'^ with gcd(a, q) = 1, we have 

J2 A{[an + f3\) = a-^ ^ A{m) + O {N^-^) 

rn^N m^laN+/3\ 
Ian+I3\=a (mod <j) rn=a (mod g) 

where the implied constant depends only on a and (3. 



Corollary 5.5. Under the conditions of Theorem 5.4, for any constant B > 
and uniformly for all integers N ^ 3 and ^ a < q ^ {logN)^ with 
gcd(a, g) = 1, we have 

J2 Ai[an + p\) = ^ + 0(^N exp (-C ^log ) ) 

lan+(3\=a (mod q) 

for some constant C > that depends only on a, (3 and B. 

Corollary 5.6. Suppose that {a,q) = (0,1) or {a,q) = (1,2). Then, under 
the conditions of Theorem \5.4\ for any constant B > and uniformly for all 
integers N ^ 3, we have 



J2 A{[an + f3\) = N + {N exp (-c(log Ar)^/^(loglogiV)~^/^)) 

[ora+/3J=a (mod q) 

for some absolute constant c > 0. 
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